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Abstract 



We re-assess constraints from flavour-changing neutral currents in the kaon system 
on supersymmetric scenarios with a light gluino, two heavy generations of squarks 
. and a lighter third generation. We compute for the first time limits in scenarios 

with three heavy squark families, taking into account QCD corrections at the next- 
to-leading order. We provide ready to use magic numbers for the computation of 
the Wilson coefficients at 2 GeV for these scenarios. 
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1 Introduction 



The initial searches for supersymmetric particles at the LHC indicate that minimal su- 
persymmetric scenarios where all supersymmetric particles have similar masses below 
a TeV are not realized in nature. This, in addition to stringent bounds from super- 
symmetric contributions to flavour- changing neutral currents (FCNC), makes scenarios 
where some supersymmetric particles are significantly heavier than others more appeal- 
ing, although these scenarios have been studied already for over a decade. With this 
motivation in mind and following the work in [I], we re-assess limits from FCNC in 
scenarios with two or three heavy families of squarks {rrig < 10 TeV), while keeping the 
gluino mass below 2 TeV. Studies of the QCD corrections in such setups, in particular 
for two heavy families, were performed in [2], [3] and [ , ]. The first work consid- 
ered leading-order (LO) QCD corrections and used the vacuum insertion approximation 
(VIA) for the hadronic matrix elements. The second one calculated next-to-LO (NLO) 
corrections and took into account lattice results for the bag parameters appearing in 
the matrix elements. More recently, [1] discussed a specific pattern of family symmetry 
breaking with hierarchical squarks and pointed out the need for a careful treatment of 
certain box diagrams in some cases, which was studied in more detail in [5]. 

The most dangerous supersymmetric contributions to FCNC often occur in the kaon 
sector, in particular contributions to Amx and the CP violation parameter e. Our 
motivations for revisiting the corresponding limits are (i) significant progress in the 
determination of the experimental values and theoretical expectations in the Standard 
Model (SM) since the publication of [3], and (ii) the need for an NLO calculation of the 
QCD corrections in the case of three heavy squark generations. 

We proceed as follows: in section 2, we summarize, for the sake of clarity, the 
way NLO QCD corrections are addressed for the AS* = 2 processes involving light 
gluinos and heavy scalars. We provide formulas that can be used to compute easily the 
renormalization group (RG) evolution of the Wilson coefficients relevant for AS* = 2 
processes, both for the case of two heavy squark generations and for the case of three 
heavy families. We compare our results with the literature, finding that in some cases 
a change of basis was missing in earlier work. In section 3, we re-assess limits from 
ArriK on the mass rriq of the heavy families of down-type squarks, for given values of the 
fiavour- violating parameters in the mass insertion (MI) approximation. We comment on 
the difference between the scenarios with three and two families of heavy scalars. In this 
section, we also determine bounds on the real and imaginary parts of fiavour-violating 
parameters for given values of rrig and the mass of the gluino, rrig. 
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2 QCD corrections for heavy squarks and a light 
gluino 

2.1 Renormalization group evolution of the Wilson coefficients 

The effective Hamiltonian for AS = 2 transitions can be written as 

5 3 



(1) 



1=1 



i=l 



where the operators Oi are 

O3 = d^'Pi^s^d^Pi^s", O4 = rPLs'^d^Pns^, 

O5 = rP^s^d^Pns'', 



(2) 



and Pl,r are chirahty projection operators, a, (3 are cofor indices, Oi = Oj (L ^ R), 
{Oi) = (Oi), and Ci = Cj (L <H- R). The RG evolution of the Wilson coefficients from 
a high-energy scale M to a lower energy /x is determined by the 5x5 evolution matrix 
IV[/i,M], 



Cifi) = W[fx,M]CiM), 



(3) 



where C represents a column vector with the 5 components Ci. We denote row vectors 
with an arrow, such as the row vector O built up by the components Oi. At LO and 
NLO, respectively, W[fi, M] can be expressed as [6] 



l^[/i,M]LO 
l?[/i,M] NLO 

/3o 



U[f^,M] 
' as{M) 

U[fi,M] + ^^J{nf)U[fi,M] - ^^^^U[fi,M]J{nf) 



2nf), 



47r 



(4) 



where Nc = 3 is the number of colors and nf equals the number of active fermion flavors, 
Hf, at energies below the gluino mass, rrig. At higher energies. 



(5) 



where riq is the number of light squarks with a mass similar to rUg [3]. The 5x5 
matrices J were calculated in [6]. They depend on the renormalization scheme. We use 
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the LRI scheme because it is also used in the lattice determination of the low-energy 
matrix elements. The one- loop anomalous dimension matrix (ADM), 7*^°), is scheme- 
independent. It depends only on Nc and is a 5 x 5 matrix, since the five operators Oi 
entering the Hamiltonian (1) do not mix with others [. ] during the evolution down to 
low energies. 

In [6], the ADMs and J are given in the Fierz basis of operators O^, defined by 



O5 



1 



ot - -Oi 



(6) 



Then we can define a matrix V to transform between the basis of the operators Oi and 
the one of the operators Of , 



O 

C 



V 



/I 





1 

\ 









1 






\ 





0/ 



(7) 



We will denote all evolution matrices in the Fierz basis with a tilde instead of a hat, 



l^[/i,M] 



v-^u[^Ji,M]v, 

y~^l?[/x,M]V. 



(8) 



In the scenario under consideration, the RG evolution of the Wilson coefficients starts 
at the mass scale of the heavy squarks, M = rriq, and ends at /i = 2 GeV, where the 
matrix elements of the AS" = 2 operators are calculated. Along the way, the gluinos and 
heavy quarks have to be integrated out at their mass scales. In the case of two heavy 
squark generations, we assume the lighter squarks to have a similar mass as the gluino. 
Hence, they are integrated out at nig, too. The corresponding mass scales are depicted 
in Fig. 1. The figure also shows the matrices governing the evolution between the scales. 
The NLO evolution from rriq down to /i yields 



C(/i) = V-'W[fi,mg]VC{mg), 
W[n,mq] = W[iJ,,mb]W[mb,mt]W[mt,mg]W[mg,mg] 



NLO 



(9) 
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:= W[mg,mg\, Ui, J(10) or J(9) 


m - 




■.= W[mt,mg], U2, J(6) 
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■.= W[m,,mt], Us, J(5) 


fi- 




■.= W[fi,m,], f/4, J(4) 



Figure 1: Energy thresholds, evolution matrices and J{nf) involved in the RG evolution of the Wilson 
coefficients. In the energy range between nig and rriq, J(10) and J(9) refer to the case of 2 and 3 heavy 
squark generations, respectively. In W4 we take fj, = 2 GeV. 



The LO evolution is analogous, with the replacement W ^ U. The subscript "NLO" 
reflects the fact that we have to truncate the expansion of the product of evolution 
matrices such that it contains only terms up to 0{as) [ ]■ Explicitly, 

W[fx, rrig] = UjjsU2Ui + UjjsU2{Wi - U^) + Ujjs{W2 - U2)Ui 

+ U^{Ws - Us)U2Ui + (1^4 - UaWsWi. (10) 

W[n,mq\ depends on the values of as at the scales rriq, rrig, nit, and /i. Plugging 
in the known values of as{mt), as{mb) and we can write the evolution matrices 

W[fi,mg] and W[mg,mq] as functions of the free parameters ^^(m^) and as{mq) as well 
as the so-called magic numbers, analogously to Eqs. (9) of [3], 
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W[mg,mq] = ^[c^M + r/eeM+r^-r^e/Mjr/f, (11) 

where 



r=l 



m = I \ , Vg = I \ , (12) 

and the magic numbers h'^'^\ c^'^\ e*^^-* and f^^^ are 5x5 matrices. However, as the 
product W[^, mg]W[mg, rriq] contains terms of order a^, to be consistent, we have to use 
the single evolution matrix 

W[^,, mq] = J2 { Ko + ^N^o + Vec^'^] d^'^ + 6^ [e^ + r^,/(^)] vAvtvt (13) 



rather than Eqs. (11). This approach is useful because for a given model, with three or 
two heavy families of squarks, we have to calculate only 1]^ and rjg. Together with the 
magic numbers, Eq. (13) then immediately yields the values of the Wilson coefficients 
at /i = 2 GeV. 



2.2 Evolution between squark and gluino mass scales 

Due to the hierarchy between gluino and squark masses, integrating out all superparticles 
at the same scale would produce large logarithms. Therefore, we proceed in two steps, 
first integrating out the heavy squark generations at rriq [2, 3] and then the gluinos and 
fight squarks at nig. 

The matching at nig is visualized in Fig. 2. In the full theory, FCNC in the neutral 
kaon system stem from the AS = 2 box diagram (I), where the scalar lines represent 
squark mass eigenstates. Using the MI approximation, we work with flavour eigenstates 
and flavour- changing mass insertions, represented by dashed lines and crosses in dia- 
grams (IIa,b), and consider only diagram (Ila). This diagram contains only d and s 
squarks. Squarks of the third generation appear in box diagrams with at least three MI, 
i.e., in diagram (lib) and higher orders. Thus, they do not contribute to the matching in 
the considered approximation. Consequently, even in the case of light third-genaration 
squarks, there are no contributions from diagrams involving one heavy and one light 
squark, which require a special treatment [1, 5]. ^ 

In the effective theory below rriq but above rrig, diagrams (Ilia) and (Illb) of Fig. 2 
yield FCNC. The former contains the A^" = 2 operators Oi and Oi, which are repre- 
sented by the dot, and it is suppressed by l/m?. The latter diagram is suppressed by 
m?/m| and thus negligible [ ]. 

Altogether, the final matching condition is (Ila) = (Ilia). The resulting expressions 
for the gluino contributions to the Wilson coefficients at high energy, Cf (m^), are given 
in Eq. (6) of [l] in terms of the MI parameters 6'^ and in the notation we use.^ They were 
calculated in [9, 2] and are the same in both cases for the squark masses we consider. 

From rrig down to rUg, gluinos and (if present) light squarks influence the RG evo- 
lution of the Wilson coefficients by changing the running of and the two-loop ADM 
of the AS = 2 operators. The impact on the ADM stems only from loops with gluinos 
and light squarks in the gluon propagator. Consequently, it can be taken into account 
by using J{nf) of [(>], where only loops with fermions and gluons were considered, and 
replacing nj by Uf [•>], as mentioned above. Explicitly, = 9 if all squarks are heavy. 
If only two squark generations are heavy, rig = 4 squarks are light and thus Uf = 10. 



^This can be different in scenarios with additional assumptions tliat cause CP violation to originate 
only from diagrams involving third-generation squarks, for example, Minimal Flavour Violation [-1]. 

^We define := 2 2 ^^^^^ X,Y e {L,R} and where one has to use the soft 

V dXY-'^^^ dXY^JJ 

mass squared matrices in the super-CKM basis, where Yukawa coupHngs are diagonaL 



5 



(I) 



(Ila) 



(lib) 







+ 



+ . . . 



-x- 

(Ila) 



(Ilia) 



(Illb) 



Figure 2: Diagrams representing the process of integrating out the heavy particles. The series of 
diagrams in the upper hue represents the transition from the complete AS — 2 box diagram to the mass 
insertion approximation, where the leading order is represented by diagram (Ila) and the next order 
by diagram (lib). Crosses denote mass insertions. Diagrams (Ilia) and (Illb) arise from diagram (Ila) 
when integrating out the d-type squarks. 

Using these considerations and the appropriate formulas from the previous subsection, 
we calculated the magic numbers determining the NLO QCD corrections in both cases. 
We list them in Appendix B. Our results differ from those of [■>]. As the disagreement 
in the magic numbers d, e and / vanishes for two heavy squark generations if we set 
V = 1, we conclude that [6] did not take into account the change of basis (7). 

After the gluino and possible light squarks decouple at rrig, the system behaves as 
in [6, 7, 8]. Therefore, the magic numbers b and c should equal those given in [7], up 
to a different ordering of the components r and small differences due to changes in the 
experimental input parameters. We find larger differences than expected for some of 
the magic numbers. The main reason seems to be that 0{a1) terms were not discarded 



3 Limits on flavour- and CP-violating parameters 
from ArriK and e 

3.1 Limits on the squark masses from Attik 

In order to set a limit on the supersymmetric contribution to ArriK, we write 



m \ i 



Ami = ^rnf + Ami 



(14) 
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where AmJ^ represents the total theoretical value of Aitlk, while Am^ and Am^j^, re- 
spectively, represent the SM and the supersymmetric contributions due to g-q diagrams. 
Even in the SM, precise computations of Am^ are not possible due to unknown long- 
distance contributions [HI]. Hence, the best one can do is to compare the best estimate 
obtained from the short-distance contributions [11], denoted by Am|P, to the experi- 
mental value 

Am^P = (3.483 ± 0.0059) x 10"^^ GeV, 

Amf = (3.1 ± 1.2) X 10"^^ GeV, (15) 

where all quoted errors correspond to the 1 a CL. Am|P was obtained by taking into 
account NNLO contributions from the charm quark. We can see that its central value 
already accounts for an 86 % of the experimental central value, but its uncertainty can 
easily account for the reported experimental value within the 1 a CL. Given the lack 
of information on long-distance contributions, the best we can do to extract limits in 
extensions of the SM is to use Am^^ — Am^ as a constraint on the order of magnitude 
of Am|^. In short, it cannot exceed 10~^^ GeV, but lower limits based on whether or not 
Am^ + Am^j^ can actually be as large as the experimental value cannot be obtained. 
Furthermore, in some cases the importance of NLO corrections cannot be really appre- 
ciated because the differences between LO and NLO lie in the range of some units of 
10~^^ GeV. Therefore, we have to be cautious when using this for setting bounds in 
models for physics beyond the SM. As in the case of e', there are many models that can 
pretty easily yield a change of Amx by some units of 10~^^ GeV, so meaningful con- 
straints can indeed be obtained. However, we have to include a generous consideration 
of all the possible uncertainties in the calculation of Am|P = 2 Re {K^\H^^'^\K ). The 
form of this theoretical expression at the known QCD level can be found in [11]. For 
the gluino contribution, we have 

5 3 

Ami = 2 Re {K''\H^^=^\lf) = 2 Re CfO, + ^ Cf 0,|¥°) . (16) 

i=l i=l 

The Wilson coefficients at /i = 2 GeV are determined by Eq. (13). The matrix elements 
of the operators are [7] 

{K'\0^\K') = ^MKflBM, 

{K'mlf) = kJ \ MKflB,{^^), ^ = 2,3,4,5, 

k^ = ^(-|>i2,i)r (17) 

These expressions differ by a factor of from the corresponding ones in lattice studies 
such as [12, 13, 14, 15]. A factor ^ stems from a different definition of the operators 
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and a factor from the normalization of the kaon states. The definition (17) ensures 
that the numerical values of the bag parameters Bi are the same in both conventions. 
The VIA corresponds to 



(K°|Oi|ir ) 



VIA 



VIA 



K 



{K^\Oa\K) 



VIA 



ms(/i) + md(/i) 
Mk 



1 1 

24 4 Vm,(/i) + md(/i) 



2 



VIA 



1 



M 



K 



12 \msilJ,) + rUdilJ,) 



i — 2,3, ki 



V 3' 3/? 



(18) 



In Appendix A we give the values of the experimental and the lattice parameters we 
use. 

With the above points in mind, once we have the value of Am^, its uncertainty, 
that we call c^msp, the experimental value and its uncertainty, a/^ 

limit on the supersymmetric contribution, Am^, using 



^cxp^ we can set a 



|Amf,| < Am^P - Amf^ + 2 (a^^sp + cta. 



2.8 X 10-^^ GeV, 



(19) 



where Am|P and Am^^ denote the central values. 
In order to compare to [ ], we consider the cases 



I : 




12 — 




12 = {sD 


12 




12 


= 0, 


II : 




12 = 




12 = (4l) 


12 


= (4r) 


12 


= 0, 


III : 


(^ll) 


12 = 


(^rr)i2 = 




12 


= 


12 


= 0, 


IV : 




)l2 = 


('^Rl)i2 




)l2 




)l2 


= 



(20) 

In Fig. 3, we present the lower limits on the heavy squark masses rriq according to Eq. (19) 
as a function of the mass of the gluino m^, for a fixed value of the flavour parameter 
K = 0.22. The plots in Fig. 3 serve mainly (i) for comparison to Figs. 1 and 2 of [3]^, 
(ii) to re-assess the importance of the NLO QCD corrections and (iii) to understand the 
difference between the scenarios with two and three heavy families of squarks. We can 
see that the latter difference is not significant. In principle, this justifies using the QCD 
corrections of [ '] also for scenarios with three heavy squark families, as done in [1] as a 
first approximation. However, as mentioned above, a change of basis was missing in [3]. 
Therefore, the limits shown in Fig. 3 are different from those in [:>]. Nevertheless, the 



^This is why we also plotted values of rrig lower than the LHC bound around 900 GeV. 
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Figure 3: Lower limits on the squark masses for the different cases of Eq. (20) with K — 0.22, derived 
from Eq. (19). The blue continuous line corresponds to the NLO evolution, the medium-dashed (red) 
line corresponds to the LO approximation, the dotted (orange) line corresponds to the LO evolution 
using the VIA for the bag parameters and the dot-dashed (light blue) line to the results without QCD 
corrections. We have indicated the labels just for case I, but for the rest of the cases the coding is the 
same. For comparison, we have additionally plotted a long-dashed black line which represents the NLO 
evolution for the case where only two squark families are heavy. 



qualitative behavior and order of magnitude of the hmits are the same. Regarding the 
QCD corrections, each step of improving the accuracy (considering LO QCD corrections 
with VIA, including lattice bag parameters and finally taking into account NLO QCD 
corrections) can have a drastic effect. Only case II may be considered an exception, 
since here some of the corrections happen to cancel partially. We can also see that NLO 
QCD corrections, as opposed to just LO corrections, are indeed relevant, especially for 
cases III and IV. 

We note that using different values of K can lower or raise considerably the limit 
on rriq. For example, for case III and K = 0.1, the limit becomes roughly rrig > 27TeV 
[16], as opposed to rriq > 56 TeV for K = 0.22. 
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Figure 4: Upper limits on y |Re(((5^y)i2((5^/y/)i2)|, derived using AniK- In all cases only the flavour- 
violating parameters shown on the vertical axis have been set to non-zero values. The limits on 
y^|Re((5pf;)^2l ^-nd ^|Re((5^L)^2l can be obtained by interchanging R o L in the upper two panels. The 

two lower panels show, respectively, the cases, where ^J\^e{{Sf^l^)l2{S^Y^)l2)\ or y^|Re(((5Lp)i2(5pL)i2)l 
are not zero. The different curves correspond, from bottom to top, to rriq — A, . . . , 10 TeV. 



3.2 Limits on mass insertion parameters 

Prom ArriK In Fig. 4 we present curves for fixed values of rrig, plotting rrig against 
the upper limits on ^\Re{{6j^Y)i2{Sx'Y 

,)i2)| determined by requiring a sufficiently small 
gluino contribution to Attlk, Eq. (f9). The different curves correspond, from bottom to 
top, to rriq = 4, . . . , 10 TeV. For rriq = 2 TeV and the maximum value of the gluino mass 
we use, rrig = 2 TeV, the condition rriq ^ rUg defining the scenario under consideration 
would be violated. 

The use of ArriK indeed yields relevant limits on all MI parameters (5xy)i2- They 
turn out to be strongest if Re{Sf^j^)i2 or Re((5p^L)i2 are non-zero. For cases when only 
the real part of one parameter is allowed to be non-zero and for gluino masses above 
900 GeV, the limit is typically of order 10~^ and never below 0.02. Of course, the 
lower TTiq, the lower the limit. For the cases where either ■\/|Re((5LL)i2('^RR)i2)| or 
Y/|Re((5LR)i2('^RL)i2)l are non-zero, the limits on the combinations are of C(10~^). 
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Prom e As it is also well-known, we can extract limits on the imaginary parts of some 
combinations of the parameters (5^y)i2 using the CP- violating parameter e. Currently, 
its experimental and SM [11] values are, respectively. 



|e| = (2.228 ± 0.011) X 10-^ 
^NNLO| _ (1.81 ± 0.28) X lO-l (21) 



Of course, in our case 



|Im (irO|i7^^=2|Z°) 



r (22) 
v2 Am/^ 

where E^^='^ = H§^='^ + H^^='^ with H^-'^='^ as given in Eq. (16). The limits are 
obtained using a condition analogous to Eq. (19), with the obvious replacements, yielding 

|e^| < 1.0 X IQ-^. (23) 

Varying the mass of the gluino for fixed values of rriq, we found the upper limits 
on iy|Im((5^y)i2((5|^/y/)i2)| that are presented in Fig. 5. As expected, the limits from 
e prove to be indeed very strict. If only the parameter y^|lni(5^^J^ is not zero, the 
limits, varying according the mass of the squarks, are of order 10~^. If only i^/|lni(5^^J^ 
is not zero, the limits decrease to O{10~^). If two MI parameters are non-zero, the 
corresponding limit is of order 10~^ and always smaller than for only one non-vanishing 
parameter. Again, the lower the mass rriq, the lower the limits. 



4 Conclusions 

Supersymmetric scenarios involving light gluinos and heavy squarks have lately regained 
a considerable amount of attention, since they may be a way to realize supersymmetry 
in agreement with the LHC exclusion limits and with an acceptable amount of tuning. 
As it is also well-known, in such scenarios dangerous FCNC can occur especially in the 
kaon sector. In an earlier work [ ], we considered a case with three heavy families of 
squarks with similar masses above 20 TeV. In order to estimate the gluino contributions 
to ArriK and the CP-violating parameter e, we used the evolution of Wilson coefficients 
from reference [o] , where a scenario with two heavy families of squarks and a lighter third 
family with a mass comparable to that of the gluino was considered. In order to quantify 
precisely the accuracy of this estimate, we have reviewed the way NLO QCD corrections 
are addressed in supersymmetric scenarios involving heavy squarks and light gluinos. 
We have calculated the renormalization group evolution of the Wilson coefficients of the 
AS* = 2 operators at NLO for scenarios with two and three heavy squark families. This 
is codified in Eq. (13). 

The SM determinations of Aitlk and e have improved significantly over the last 
decade. Together with important improvements on the lattice QCD bag parameters 
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Figure 5: Upper limits on ^ \lm{{SxY)i'^i^'x'Y')i^)\- The limits on y'|Im((5ppj)^2l ^'i^d ^^1^11(5^^)121 
can be obtained by interchanging R -H- L in the corresponding upper two panels. The lower two 

panels correspond to the cases where only the combinations, respectively, \/|Ini(((5LL)i2((5{^j^)i2)| or 



|Ini(((5Lf{)i2(5RL)i2)| are non-zero. The different curves correspond, from bottom to top, to 
. . . , 10 TeV. 



entering the effective Hamiltonian of tlie AS* = 2 transitions, tliis proves relevant for 
re-assessing tlie limits on supersymmetric scenarios contributing to Amx and e. This 
was the second motivation for this work. 

We have then determined the lower limits on the mass of the heavy squarks, rriq, 
coming from Amx, using the MI approximation with a value of ^/ \^g{{^xy)i2)\ = 0.22, 
X, y G {L, R}, for one or two flavour- violating parameters. For cases with only one non- 
zero flavour- violating parameter, the limits on rrig for a gluino mass of about 900 GeV 
and three heavy squark families are as follows: if only either -\/|Re((5^^jf^ = 0.22 
or A/|Re((5^j^)^2) I = 0.22 is non-zero, rriq has to be larger than 5 TeV. If only either 
V|Re((5j^R)?2)| = 0.22 or v^Re((5^jy] = 0.22 is non-zero, the limit is above 14 TeV. 
If only the first two squark families are heavy, the limit changes by about 100 GeV. For 
the case of v^MI^SSM^)^ = 0-22 and the case of ^/\M(^tJ^MJ^\ = 0-22, 
while the other fiavour-violating parameters are zero, the limits are respectively around 
56 TeV and 125 TeV, again for gluinos around 900 GeV. For two heavy families of 
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squarks, the limits are a few TeV larger than for three heavy families. We conclude that 
when the gluino contribution to Aitlk is close to its upper limit, it changes only by a 
few percent if the mass of the third-generation squarks varies between rrig and rriq. 

We have also noticed that indeed the NLO QCD evolution, as opposed to just the 
LO evolution, proves relevant, especially for cases where there is more than one type of 
flavour- violating parameter (Sxy) 12 involved [ ]. As we have seen, the difference can be 
as large as 20 %. 

Finally, we have obtained bounds on the real and imaginary parts of the combinations 
of flavour- violating parameters (5^y)i2 that enter into the Wilson coefficients Cf and Cf 
for three heavy families of squarks. The former bounds stem from Attik and the latter 
from e. We have explored these limits for squark masses between 4 and 10 TeV, varying 
the mass of the gluino for each value. In short, when only the real part of one type of 
parameter (Sxy) 12 is allowed to exist and for masses of the gluino above 900 GeV, the 
limit on |Re(5^y)i2| is at most of O{10~^) in most cases, of course the lower the mass 
rriq, the lower the limit. When only ^|Im(5^y)^2l type of parameter is allowed 

to exist, its limit is of O{10~^) or smaller. Of course, when several types of flavour- and 
CP- violating parameters (5xy)i2 are non- vanishing, the upper limits become smaller. 
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A Experimental information 

Values of experimental and lattice QCD parameters not given in the main part are 
listed in Tab. 1. We are aware that there have recently been some improvements in the 
determination of [17, 18, 19, 20, 21, 22], which have been averaged to a value of 
0.7643 ± 0.0097 [23]. This would change the values for Amx and e obtained in [ii], but 
we have checked that the impact on our results is negligible, since the limits are mainly 
determined by 0"AmA' and a^- 
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Experimental & lattice inputs 





0.1184 ± 0.0007 


[ ] 


mt := Mf°'^ 


(173.5 ± 0.6 ± 0.8) GeV 


[ ] 


rUb := mbirrib) 


(4.18 ± 0.03) GeV 


[ ] 


Ik 


(0.1561 ± 0.00085) GeV 


r "1 

[ ] 


Mk 


(0.497614 ± 0.000022) GeV 


[--] 


Bi 


0.51 ±0.02 


[15] 


B2 


0.73 ± 0.04 


[15] 


53 


1.29 ±0.11 


[15] 


54 


1.04 ±0.07 


[15] 


B5 


0.76 ±0.09 


[15] 


m,(2 GeV)^i 


(0.12 ±0.006) GeV 


* 


mrf(2 GeV)^^ 


(0.006 ± 0.001) GeV 


* 



Table 1: Experimental and lattice QCD values used for our analysis. The bag parameters Bi are given 
in the RI scheme. *We have computed with RunDec [2'i] the values of and irid in the RI scheme, 
using as input the values given by the PDG [ -1] in the MS scheme. 

B Magic numbers 
B.l Evolution below rrin 



8 1 1 - V24T 1 ± V24T 
'7' 7' 21 ' 21~ 



(0.29, -1.1, 0.14, -0.69, 0.79)^ (24) 



"loii 

"L023 

dr) 
'^L032 

dr) 
'^L033 

"L044 

'^L045 

"l055 

dr) 

"NLOll 
dr) 

"NL022 
"NL023 
"NL032 



(0.77,0,0,0,0)^ 

(0,0,0,1.8,0.0083)^ 

(0,0,0,-0.48,0.13)^ 

(0,0,0,-0.12,0.032)^ 

(0,0,0,0.032,0.48)^ 

(0,2.8,0,0,0)^ 

(0,0.94,-0.29,0,0)^ 

(0,0, 0.88, 0,0)^ 

(0.045,0,0,0,0)^ 

(0,0,0,0.51,0.0020)^ 

(0,0,0,-0.13,0.030)^ 

(0,0,0,0.088,-0.0036), 



(25) 
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'^NLOSS 

"NL044 
llr) 

"NL045 
"NL054 



iSr) 

"NL055 
^11 



C 



(r) _ 

22 ~ 

Jr) _ 

'-23 — 

Sr) _ 

^32 ~ 

Ir) _ 

^33 ~" 

(r) _ 

C44 — 

Sr) _ 

'^45 — 

^(0 _ 

<-54 — 

^(0 _ 



(0,0,0,-0.023,-0.055), 

(0, 1.3,0,0,0),. 

(0,0.42,0.096,0,0),. 

(0,0.14,0,0,0),. 

(0,0.048,-0.063,0,0),. 

(-0.015,0,0,0,0),. 

(0,0,0,-0.19,-0.0026), 

(0,0,0,-0.015,0.0063),. 

(0,0,0,0.012,-0.0099)^ 

(0, 0, 0, 0.00096, 0.024)^ 

(0,-0.51,0.0054,0,0),. 

(0,-0.26,-0.0061,0,0), 

(0,0,-0.016,0,0),. 

(0,0,0.018,0,0)^ 



(26) 



(27) 



B.2 Evolution above m? 



"44 

"45 



d 



(r) 
55 



(1.0,0,0,0,0),. 
(0,0, 0,0.98,0.017)^ 
(0,0,0, -0.26, 0.26)r 
(0,0,0,-0.064,0.064), 
(0,0,0,0.017,0.98),. 
(0,1.0,0,0,0),. 
(0,0.33,-0.33,0,0),. 
(0,0, 1.0,0,0),. 



(28) 



These magic numbers correspond to LO effects. Consequently, they are independent of 
77,/, in particular of the number of heavy squarks. 



B.2.1 Three heavy squark generations 



1 

5' 



1 - ^241 1 + \/24T^ 



15 



15 



(0.40, -1.6, 0.20, -0.97, 1.1),. (29) 



15 



^11 




(0.0019,0,0,0,0)^ 






(0,0,0,0.15,0.00044)^ 


Jr) 
^23 




(0, 0, 0, -0.040, 0.0067)^ 


Jr) 
^32 




(0,0,0,0.013, -0.0059)^ 


^33 




(0,0,0,-0.0034,-0.091)^ 


Jr) 


— 


(0,0.27,0,0,0)^ 


Jr) 
^45 


— 


(0,0.091,0.037,0,0)^ 


Jr) 




(0,0.018,0,0,0)^ 


^55 




(0,0.0061, -0.035, 0,0)r 


Ar) 
Jll 


= 


(-0.0019, 0,0,0,0)^ 


Ar) 
J 22 




(0,0,0,-0.15,-0.0044)^ 


Ar) 
J 23 




(0,0,0,0.0088,0.025)^ 


Ar) 
J 32 




(0,0,0,0.0098,-0.017)^ 


Ar) 
J 33 




(0,0,0,-0.00058,0.095)^ 


Ar) 
J 44 




(0,-0.28,0.0061,0,0)^ 


Ar) 
J 45 




(0,-0.12,-0.0098,0,0)^ 


Ar) 
J 54 




(0,0,-0.018,0,0)^ 


Ar) 
J 55 




(0,0,0.029,0,0)^ 



(30) 



(31) 



B.2.2 Two heavy squark generations 



24 3 1 - v/24T 1 + ^241 
13' 13' 13 ' 13~ 



(0.46, -1.8, 0.23, -1.1, 1.3)^ (32) 



Jr) _ 

^11 — 

Ar) _ 

^22 — 

Ar) _ 

^23 — 

Ar) _ 

^32 — 

Ar) _ 

^33 — 

Ar) _ 

C44 — 



^r) 
--45 



(-0.012,0,0,0,0)^ 
(0,0,0, 0.19, -0.00026)^ 
(0,0,0,-0.050,-0.0040), 
(0,0,0, 0.010, -0.0084)^ 
(0,0,0, -0.0027, -0.13)^ 
(0,0.35,0,0,0)^ 
(0,0.12,0.038,0,0)^ 
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(r) 
^54 

(r) 
^55 

ill 

i23 

Ar) 
J 32 

Ar) 
J 33 

Ar) 

Ju 

Ar) 
J 45 

Ar) 

J 54: 

Ar) 
J 55 



0,0.018,0,0,0)^ 

0,0.0060,-0.042,0,0), 

0.012,0,0,0,0)^ 

0,0,0,-0.19,-0.0038) 

0,0,0,0.020,0.035)^ 

0,0,0,0.012,-0.014)^ 

0,0,0,-0.0013,0.13)^ 

0,-0.35,0.0060,0,0)^ 

0,-0.14,-0.012,0,0)^ 

0,0,-0.018,0,0)^ 

0,0,0.036,0,0)^ 



(33) 



(34) 
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